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Abstract 

Few facts are known about the entanglement entropy for disconnected regions in quantum field 
theory. We study here the property of extensivity of the mutual information, which holds for free 
Q\ . massless fermions in two dimensions. We uncover the structure of the entropy function in the extensive 

case, and find an interesting connection with the renormalization group irreversibility. The solution is 
a function on space-time regions which complies with all the known requirements a relativistic entropy 
function has to satisfy. We show that the holographic ansatz of Ryu and Takayanagi, the free scalar 
Oh. and Dirac fields in dimensions greater than two, and the massive free fields in two dimensions all fail 

£^ j to be exactly extensive, disproving recent conjectures. 

1 Introduction 

The entanglement entropy enclosed in a region V of the space is ambiguous in quantum field theory 
(QFT). This is because there is a divergent amount of entanglement between the V and the exterior, 
due to the vacuum fluctuations. However, as it happens with the Casimir effect, the unwanted effects of 
^ . these fluctuations can be governed by looking at different, cutoff independent quantities, which can be 
constructed from the local entropy. With this aim let us introduce the mutual information between two 
CN ■ non intersecting regions A and B, defined by 
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I(A, B) = S(A) + S(B) - S(A U B) , (1) 



> 

m 
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where the entanglement entropy S(V) = —trpv log py is the the von Neumann entropy of the local density 
matrix py corresponding to the global vacuum state. Since the divergences in S(V) are related to the 
boundary of V, I (A, B) is free from divergences and ambiguities, and thus a good measure of entanglement 
in QFT. 

The calculation of I(A, B) requires the knowledge of the entanglement entropy for multicomponent 
sets, from which very little is known so far. In particular, the only exact result is for free massless fermions 
in two dimensions. For a Dirac field the multicomponent entropy function is [UH], 



S(X 1 U X 2 U ... U X p ) = - j log |oi - bj\ - log h - aj\ - log \bi - bj\ - p log e J , (2) 

where ctj and bi are the left and right endpoints of the interval Xi and e is a short distance ultraviolet 
cutoff. Remarkably, this specific form of the entropy gives an extensive mutual information [2] , 

I{A,BUC) = I(A,B)+I(A,C). (3) 

Given this peculiar result, one may naturally wonder if there are other QFT with the same property. 
For a couple of years it was thought that the expression ^§ for the entropy of the free massless Dirac field, 
multiplied by a global factor of the Virasoro central charge, was valid for any two dimensional conformal 
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field theory (CFT) [3j. Recently, the derivation of this result was called into question (see [I] and the 
note added to [3]), and up to the present, the expression of the entropy for more than one interval and 
more general CFT remains unknown. The work |5] reports numerical results for the compactified scalar 
which hint to a different behavior than ([2]). 

In this work, we will investigate the property of extensivity of the mutual information in QFT from 
a general point of view. The results allow us to disprove extensivity for some theories in more than two 
dimensions. In the course of the investigation we also find an unexpected connection with renormalization 
group irreversibility. 

Two recent developments add interest to this subject. The first one is the geometric ansatz for the 
entanglement entropy of a CFT proposed by Ryu and Takayanagi in the context of the Maldacena duality 
[BJ. This gives the entanglement entropy of a connected spatial region V on the Minkowskian boundary 
of an AdS space as the area of a minimal surface in the bulk, whose boundary coincides with the one of 
V. This holographic entropy has been supported by several arguments [TUB]. However, this construction 
does not provide a formula for the entropy corresponding to multicomponent sets. In [9] Hubeny and 
Rangamani propose one such extension (see also [6j). For the two dimensional case it is a generalization 
of © 

s(x 1 ux 2 u ...ux p ) = J2Si(\ai - bj\) - J>i(k - ojl) - Y, Sl ^ - b iD ' ( 4 ) 

i,j i<j i<j 

where S\(l) is the entropy function for single component sets. This formula leads also to extensivity of 
the mutual information. For more dimensions the proposal in [9j gives an extensive I (A, B) at least for 
a particular relative position of the sets involved. 

A different kind of investigation has also met with the present problem |10j . There, we have found 
that an unambiguous definition of the entropy contained in the Hawking radiation on a finite region V 
is given by half the mutual information between the black hole and V. We have also argued using the 
mutual information that, in case there exists a Hagedorn transition, there is a natural regularization of 
the entanglement entropy in QFT which does not requires quantum gravity. This may hint to a different 
understanding of the finitness of the black hole entropy in terms of the maximum ammount of information 
that can be shared between the inside and the outside of the horizon. The interpretation of the ordinary 
thermodynamical radiation entropy (the entropy in the Hawking radiation) as mutual information requires 
the extensivity of this later at least in the particular limit case when V is large and far from the black 
hole. The conditions for this to hold are unknown at present, and it also motivates further investigations 
on the properties of the mutual information in QFT. 

2 General considerations 

The fact that I (A, B) is a well defined quantity, and that we take the vacuum as global state, has the 
added advantage that the spatial symmetries of the theory (Poincare, conformal) are explicitly manifest. 

Further relations come from causality [11]. Indeed, heuristically one would expect that the evolution 
between the Cauchy slices A and A' in figure 1 is given by a unitary operator in quantum field theory. 
This would imply that the entropy is the same for these surfaces, S(A) = S(A'). In other terms, the 
causal evolution does not allow the information present in A' to move at a velocity greater than the 
light speed and then not to be present also in A (escaping from the diamond shaped set A in figure 
1). The information content must then remain constant between these sets. There is a subtle point 
however, since the entanglement entropy requires regularization and this may spoil the unitarity of the 
evolution. However, this problem is only apparent since it is solved (at the expense of dealing with two sets 
rather than one) by the use of the mutual information, which is regularization independent. The mutual 
information can be defined (at the full rigor of the axiomatic level in QFT) using only the global state 
and the algebras of the operators pertaining to two different regions A and B. These later are exactly the 
same as the algebras of local operators corresponding to A' and B' due to the unitary causal evolution of 
the Heisemberg operators. Thus, these algebras, and the mutual information, are in fact funcions of the 
"diamond shaped" or causally complete sets |12j . These D dimensional sets, where D is the space-time 
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Figure 1: Two spatial sets A and B in two dimensional Minkowski space-time. Light rays are plotted at 
±45°. The spatial sets A, A' , and B, B', have the same causal domain of dependence A and B respectively. 
Then we have I(A, B) = I (A', B'). 

dimension, are the domain of dependence of spatial sets and in this sense they are classes of equivalence of 
Cauchy slices which have the same domain of dependence. The causality is then explicitly manifest in the 
mutual information leading to I(A, B) = I(A' , B') whenever A and A', and B and B' are Cauchy surfaces 
for the same D dimensional sets A and B. We note that A and B have to be at a space-like distance in 
order to define a mutual information, since the respective operator algebras have to be commuting. 

One can naturally wonder if all the information about universal (regularization independent) quantities 
contained in the entropy are already present in the function I (A, B). Though it is not a theorem, a positive 
answer may be expected. This is because, due to the purity of the vacuum state, we have S(A) = S(— A), 
where —A is the set complementary to A, and the entropy corresponding to the whole space A U —A is 
zero (again due to the purity of the vacuum |0)). Thus, we recover the entropy function from the mutual 
information by taking the limit 

B lun A I(A,B) = 2S(A). (5) 

In this limit the separation between A and B provides a cutoff, and the mutual information diverges when 
A and B come into contact, corresponding to a divergent S(A). 

Some properties of I(A, B) can be deduced from its relation (P) with the entropy. The strong subad- 
ditivity of the entropy [13] 

S(A) + S(B) > S{A U B) + S{A n B) , (6) 
implies that the mutual information is positive and monotonically increasing 

I(A,B) > 0, (7) 
I(A,B) < I(A,C), BCC. (8) 

In order to measure the degree of extensivity of the mutual information we can introduce the tripartite 
information function (somehow information shared between B and C with respect to A) 

I(A, B, C) = I{A, B) + I(A, C) - I(A, BUC). (9) 

This function has complete permutation symmetry since 

I(A, B, C) = S{A) + S(B) + 5(C) - S(A U B) - S(A U C) - S(B U C) + S(A U B U C) . (10) 

For extensive functions (i.e. ([2])) we have 

I{A,B,C) = 0. (11) 

One particular case when I(A, B, C) = is when the joint state for the three systems, paubuc is pure, a 
situation which does not hold in QFT unless that A U B U C is the whole space. The general condition 
which makes the mutual information extensive for a tripartite system, I(A,B,C) = 0, are unknown at 
present. 
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A particular limit of the tripartite information function is known to provide an order parameter 
for massive theories with topological order [J3]. In this case I(A,B,C) has been called topological 
entanglement entropy. In three dimensions this limit of I(A, B, C) is always negative or zero, a fact that 
follows from the strong subadditive property. 

One may be tempted to speculate that a natural generalization of the strong subadditive inequality ([6]) 
for tripartite systems is that I(A,B,C) has a definite sign. Unfortunately, this does not hold in general, 
and both signs are allowed [15]. Below we show some examples of this fact in QFT. 



3 Extensivity and the renormalization group irreversibility 

We now investigate the properties of theories which are extensive, that is, theories such that the mutual 
information satisfies the property ([3]), or equivalently I(A,B,C) = 0. If the mutual information is 
extensive, I(A, B) can be decomposed in a sum of the mutual information of infinitesimal sets covering 
the surfaces A and B. Moreover, the result must be the same independently of the way in which this 
decomposition into small sets is made. Thus, modulo details from meassure theory, we can write the 
mutual information in terms of an integral 

I(A,B)= / da x I da y j(x,y,rj AyX ,r] B) y) , (12) 
J A Jb 

where t\a and r/B are the normal vectors to the surfaces A and B. Note that j(x,y,rjA, x ,rjB t y) in this 
expression depends on the points x, y, which can be thought as the location of the infinitesimal spatial sets, 
and their orientation in spacetime, given by the normal vectors ija,x and f)B,y at these points. However, 
j(x, y, r]A, x ,i]B, y ) in equation (12) does not depend on other details of A and B since the same infinitesimal 
sets can be used to cover (and compute the mutual information of) any other pair of sets which passes 
through the same points x and y with the same normal orientations t]a,x and f]B,y 

Causality requires that I(A, B) must be independent of the Cauchy surfaces A and B for the causal 
completion A and B (the diamond shaped sets in figure 1). Now, one can make an infinitesimal transfor- 
mation of the Cauchy surfaces in such a way that the points x and y are varied infinitesimally, but the 
normal vectors t]a, x and r]B,y get macroscopically varied. In order that this type of ripples in the surfaces 
do not change the result of the integral in (fl~2j) the function j(x,y,r]A, x ,ilB,y) must be linear in r\A and 
t]b- Equation (fT2"j) then becomes 

I(A,B)= / da x rf Ax \ da y rf B J^ u (x,y) , (13) 
J A ' Jb 

for certain symmetric tensor J^(x,y) which is a function of the points x and y only. 

Now, insisting in that the result of (|13p must not depend on arbitrary variations of the surfaces A and 
B (while keeping them spatial and their boundary fixed), we can use the Gauss theorem in order to prove 
that this bi-current must be conserved, 

d tM J^(x,y)=0. (14) 

The Poincare invariance then gives 

{x-yY{x -y) v n ^. gT 



J ^ V) = V (J-,)^ G{^W) - (x _^ )2(D _ 1} F(Vi^W) • (15) 

We take the metric signature such that (x — y) 2 > for a pair of space-like separated points. The functions 
F(l) and G(l) are dimensionless. The conservation of J^ v {x — y) leads to 

(G-F)' = -(Z?-l) (2i? ~ G) . (16) 



In particular, for an extensive conformal theory F and G must be constant, what according to (|16p means 
G = 2F, in any dimension. 
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Another constraint to F and G follows from the positivity of the mutual information, which in the 
present extensive case also entails monotonicity. This implies 

^J^(x,y)^>0, (17) 
for any vectors 77 1 and rf in the positive light cone. This is equivalent to 

2F>G>0. (18) 
Then, according to (|16p . the dimensionless function C(Z), 

C(l) = G{l)-F{l) (19) 

is decreasing 

C"(i)<0. (20) 

This means there is a monotonous behavior along the renormalization group trajectories for extensive 
theories in any dimensions, a result resembling the C theorem in two dimensions [16J. At the fixed points 
it is C = F > 0. As we discuss in Section 4, in four dimensions, this quantity at the conformal point is 
proportional to the type a conformal anomaly Then the C'(l) < is consistent with the proposal that 
this anomaly plays the role of the Zamolodchikov's c- function in four dimensions [17J. 

The existence of an infrared fixed point then means that G cannot get to be smaller than F, 

C(l) > (21) 

for any I, and lim;_ +oc (F(/) — C(Z)) = 0. From (|16|) we have 

F = -(D -l^lC + C , (22) 
G = -{D- 1)- 1 IC' + 2C. (23) 

Thus, taking into account these equations, all the constraints are summarized by C{1) being positive and 
decreasing. 
Writing 

C{1) = (D- l)l 2D ^H'{l) (24) 

for some function H(l), we have 

•WO = -0A H{1) + g>lu d a 8 a H(l) . (25) 

The conservation of the current is transparent in this expression. With this we can integrate eq. (|13p to 
obtain 

I(A, B) = - f daf [ da% H(\x - y\) , (26) 
JdA JdB 

where the integration is now over the (D — 2) dimensional boundaries of the spatial regions A and B, taken 
with the same orientation. In particular, when A and B are contained in the same spatial hyperplane it 
is 

I(A,B)= I da x [ da y H(\x-y\) (ftf.ffi), (27) 
JdA JdB 

with ff A and fj B the outward pointing unit vectors normal to the surfaces OA and dB respectively. 
These boundaries depend on the causally complete regions A and B, and are explicitly independent of 
the particular Cauchy surfaces. The arbitrary additive integration constant in H is irrelevant since the 
surfaces dA and dB are closed. 

In the case of a conformal theory, the function H (I) has the particular form 

ff(0 = "jS?j, (28) 

where Co = C(0)/(2(L> — 1)(D — 2)), with C(0) the constant value of the function C{1) at the conformal 
point. In this case it can be checked that (I13p is invariant under the full conformal group in D dimensions. 
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Figure 2: In order to obtain the entropy corresponding to A, we can evaluate the mutual information 
between A and — A e , which is the set formed by the points separated from A more than a short distance 
cutoff e. 

4 Testing extensivity in different theories 

In this Section we compare the general formula (|27|) for the extensive case with some known exact 
results for specific theories. These are available in few cases, and correspond to free scalars and fermions 
[181 [TT?1 120] . or to the CFT in the Ryu and Takayanagi holographic ansatz. 

In order to make this comparison we take advantage of the property ©, which expresses the entan- 
glement entropy as a particular limit of the mutual information, and holds for pure states. In this Section 
we consider sets included in a single spatial hyperplane, and thus we can just impose a distance cutoff e 

S(A) = ±I(A,-A e ), (29) 

where — A e is the set of points on the hyperplane of A which are at distance longer than e from A (see 
figure 2). Expanding for e — > we get the entropy with its divergent and finite terms. 

4.1 Two dimensional case 

For the extensive case in two dimensions we derive the entropy from (|29[) using the mutual information 
(|27p . For a single interval we obtain 

H(l) = Sx(l). (30) 

That is, H(l) is just the single interval entropy (plus an arbitrary additive constant). Then we can write 
C(l) = IS[(1). This indicates that the condition C'{1) < coincides with the entanglement entropy 
c-theorem for this particular case [U EI] • 

Using (|30p in (|27p we obtain the expression for the multicomponent entropy, which coincides with 
the Hubeny-Rangamani proposal (jlj). Thus, this is equivalent to extensivity in two dimensions. As we 
mentioned, this formula holds for the free massless fermion. In contrast, it does not hold for free massive 
fields. Some numerical results about this fact were reported in [lj |2]. The figure 3 shows two particular 
examples. These are numerical calculations in the lattice of I(A,B,C) for a Dirac field with different 
masses. The involved sets are an interval A of size a separated by a distance d from the two adjacent 
intervals B of length b and C of length c. For this configuration I(A, B, C) only involves the mutual 
information function for two single component sets. Typically, the plotted points correspond to sizes from 
200 to 2000 lattice points. The examples shown are for a = b = c, d/a = 3, and a = d, b = c = 2/3a. The 
picture shows that the negative sign predominates, but actually both signs of I (A, B, C) are attained by 
the massive fermion. 

We have also done analytical calculations for free massive fermions. These are based on a direct 
diagonalization of the density matrix for the multicomponent massless case, and perturbations for small 
mass. The details of the derivation of these analytical results for the deviation from extensivity will be 
reported elsewhere. Here we just present the result for the mutual information of two intervals A and B 
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Figure 3: The tripartite information function I(A, B, C) for a two dimensional Dirac field with different 
masses. The involved sets are an interval A of size o separated by a distance d from the two adjacent 
intervals B of length b and C of length c. The examples shown are for a = b = c, d/a = 3 (the curve with 
the smaller well), and a = d, b = c = 2/3a (the curve at the bottom, with the larger well). The x axis is 
the product of the mass with the total length a + d + b + c. 



of lengths a and b separated by a distance d> a,b, in the small mass limit, md <C 1. We have for a Dirac 
field 

+i ((2 \og(d) + 2 lE -2 log(2) + 1) a 6 + a& (° + 6 ) + 0(cT 2 )) m 2 log(m) + 0(m 2 log°(m)) 

The first term is the conformal extensive contribution. The second one is the leading correction, which 
is still extensive. The last term contains the first non-extensive correction. In this case it gives negative 
tripartite information 

I(A,B,C) ~ -| ^m 2 log(m), (32) 

where as above a is the size of A, and 6 and c are the ones of the consecutive intervals B and C. 

For the free scalar the deviations from extensivity are much larger than for the fermion, and there 
is no extensivity in the massless limit either. This is because the homogeneous component of the field 
is a zero mode in this limit. The correlation function diverges logarithmically with the mass, and thus 
the typical size of the fluctuations on the homogeneous mode grows as (— log(m)) 1 / 2 . Correspondingly, 
the entropy grows as the logarithm of this volume in field space [22], and becomes infrared divergent 
S(A) ~ 1/2 log(— log(m)) |18j . This term in the entanglement entropy is independent of the number of 
components of A due to its infrared origin. In consequence, I(A,B) ~ I log(— log(m)) is also infrared 
divergent, what impedes the formula ([3]) to hold in the small mass limit. We have positive tripartite 
information 

I(A,B,C) ~±log(-log(m)) (33) 

in this case. In ref. [9] the equation was proposed to hold for all two dimensional QFT. As we see, 
this cannot be the case. 



4.2 Vertices in three dimensions 

The entanglement entropy S(V) associated to a polygonal set V in three dimensions has a term which 
runs logarithmically with the cutoff e 

slog(eA), (34) 
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Figure 4: Comparison between the coefficients of the logarithmic term in three dimensions for different 
theories. The picture shows, from top to bottom, the functions ss/se, sd/se and sh/se- The differences 
increase toward the origin where these ratios are ss/se = 1-078, sd/se = -980 and sh/se = -956 (not 
shown). The functions sh and se are normalized such that the coefficient of (x — ir) 2 in the Taylor series 
expansion around x = tt coincide with the ones of sd and ss- 

where A is a parameter with dimension of energy, depending on V and on the details of the theory. The 
dimensionless universal coefficient s, being an extensive and local function on the boundary, is given by 



where the sum is over the vertices Vi and where Xi is the vertex angle. The function s{x) was calculated 
for a free scalar field s = ss [IS], and for a Dirac one s = sp [20J. The calculations are based in the 
relation between the entropy and a partition function in D dimensions, where the field has particular 
boundary conditions imposed on the D — 1 dimensional set V [23]. 

For the holographic ansatz, a similar term in three dimensions was found by Hirata and Takayanagi 
[8]. We call their angular function sh{x). 

On the other hand, from (|27|) and (|29p we can calculate the regularized entropy of a plane angular 
sector and extract the logarithmic contribution to this entropy, obtaining the coefficient s = se of the 
"extensive" function. This is 



with C3 a constant corresponding to the ultraviolet fixed point. 

In order to compare the curves of the different angular functions s(x), we normalize sh and se to have 
the same quadratic coefficient as ss and so in the Taylor series expansion for x around ir. The quadratic 
coefficients for these later functions coincide. For more details see [20J. The functions sp, ss, se and sh, 
divided by se(%), are plotted in figure 4 for a range x £ [tt/2, tt]. Though they are all remarkably similar 
to each other, they differ by around 1% in this range. Therefore, we conclude none of these theories is 
extensive. Moreover, we are able to check that they are linearly independent, using the coefficients of the 
series expansion around x = tt reported in [20J. This means there is no non-interacting combination of 
them which has extensive mutual information. 

4.3 Dimensional reduction 

The non-extensive character of the mutual information for free massive fields in two and three dimensions 
is transferred to higher dimensions. In order to see this consider sets Ae 1} Be >1 and Cd 1 in space-time 
dimension D\, which are products of D2 dimensional sets Ad 2 , Bd 2 and Cd 2 times a large cube of 




(35) 



s E = 2 C 3 (1 + (vr - x)cot(x)) 



(36) 
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size L in the other D% — D2 remaining dimensions. Decomposing the D\ dimensional fields in modes 
with momentum parallel and transversal to the D2 dimension selected, the squared momentum in the 
transversal direction is interpreted as a squared mass for the D2 dimensional theory. We arrive in this 
way at a dimensionally reduced formula in the large L limit, which relates the higher dimensional function 
with the lower dimensional one, summed over different masses [181 



I Dl (A Dl , B Dl , C Dl , m) ~ qD (27r )Di-D 2 / dz z 01 ' 02 - 1 Id, (A D2 , B D2 , C D2 , Vm 2 + z 2 ) . 

(37) 

Here Y(S^) = (n + l)7r( n+1 )/ 2 /T((n + 3)/2) denotes the volume of the n dimensional sphere and q n is 
the multiplicity of the spin degree of freedom in n dimensions. In particular, when D2 = D\ — 2 = D we 
have 

d q£> 2 L 2 

— I D+2 (A D+ 2, B D+2 , C D+2 ,m) = -m ——I D (A D , B D , C D , m) . (38) 

dm qr> {2tt) 

This shows explicitly that for a free scalar or Dirac field in any dimension there is a mass such that 
the theory is non-extensive. Then, there is non extensivity for any nonzero mass, since by dimensional 
analysis I(A, B,C,m) = I(XA, XB, XC,m/X). For the massless case in more than three dimensions we 
can use directly (|37p to relate the tripartite information with the massive lower dimensional ones. The 
numerical calculations of figure 3 for the fermion fields show that these integrals are not zero in general, 
and thus the massless fermions are non-extensive in any dimension D > 2. We will show later the same 
holds for massless scalar s. 

Note that (|3T[) implies that Id(A, B,C,m) is not bounded for D > 2, even if the two dimensional 
function were bounded. 



4.4 Spheres in D dimensions 

Now we consider a conformal theory in D space-time dimensions, and focus attention on the mutual 
information between a ball bounded by a sphere A of radius Ra and the exterior B of a concentric sphere 
of radius Rb, with Rb > Ra- Let us introduce the parameter x = Ra/Rb with x £ (0, 1). 
The mutual information for the extensive case reads 

d-2 d d-2 (Ma-SIb) 



I(A,B) = C D j dn A J dn B RA~ 2 R 



PaRa-VbRbI 20 - 4 



The result for D odd is 



I(A B) - XZLlzlg f40 , 
I(A,B)- ^-3 ; (40) 



with 7 = C D \{S D - 2 )\{S D -' i ). For D even we have 



D , . , , D 



^ ' B) = ~2P^T lQ g (j+Z ) + 2C -4 (jD _ 3)! , (1 _ x 2 )D -2 . (41) 

where Pd{%) is a polynomial of degree {2D — 6). We have not been able to obtain a closed expression of 
Pd{x) valid for any dimensions. Table 1 shows Pd{x) for D < 10. 

Now we can compare these functions with the mutual information corresponding to the Ryu-Takayanagi 
ansatz. For the present case of two concentric spheres this is given in terms of an integral over the so- 
lution of a non linear differential equation [8]. The important point we want to stress here is that the 
holographic ansatz displays a "phase transition" (Gross-Ooguri transition [23] in the context of Wilson 
loop calculations), a non differentiable change in the curve I(A,B) as a function of x. This was studied 
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D 


Pd{x) 


4 


x 2 + 1 


6 


3x 6 - llx 4 - llx 2 + 3 


8 


15x 10 - 85x 8 + 198x 6 + 198x 4 - 85x 2 + 15 


10 


105x 14 - 805x 12 + 2681x 10 - 5053x 8 - 5053x 6 + 2681x 4 - 805x 2 + 105 



Table 1: the polynomial Pd(x) for even dimension D < 10. 



in detail for D = 3 and D = 4 in [8] (see also [25]), but is a general feature for any dimension. This is 
because the entanglement entropy of an annulus bounded by the two spheres is computed as a minimal 
area in AdS having the spheres as a boundary. When the spheres are near to each other the topology 
of this minimal surface is similar to a half torus, directly connecting the surfaces. In contrast, in the 
case the smaller sphere tends to collapse to a point, the surface of minimal area is necessarily one having 
two connected components, each one bounded by a different sphere. Then, at some x £ (0, 1) there is 
a change in the topology of the minimal surface, and a non analytical transition. This non analiticity 
is not present for the functions (|40p and (|4ip corresponding to an extensive mutual information in odd 
and even dimensions respectively. Thus, we conclude that, for any dimension, the entropy function of the 
Ryu-Takayanagi ansatz cannot correspond to an extensive theory. 



4.5 Universal coefficients 

The presence of a logarithmic divergent term in the geometric entropy for smooth sets in even dimensions 
is a general phenomenon. This is a consequence of an anomaly which is present for conical manifolds 
with curved singularity surface |26j. In agreement with this fact, the eq. (|41|) for the extensive case in 
even dimensions, leads to a logarithmically divergent term in the entropy for a sphere. The coefficient of 
log(e), with e the cutoff, reads 

"E = ^fe 2 • (42) 
For the holographic ansatz the coefficient of the logarithmic term for a sphere in even dimensions is [6] 

SH ~2G D N +l T[2=±)(D-2)\\ ' ( 3) 

It is interesting to note that this coefficient, which in the extensive conformal case is proportional to 
the function C, has been related to the central charges of the theory (coefficients of the energy momentum 
tensor anomaly) |6j. For a sphere in four dimensions the coefficient is proportional to the type a conformal 
anomaly [27J. 

Another interesting universal coefficient can be obtained from the spherical corona in the limit of small 
separation x = R1/R2 ~ 1. This is locally equivalent to the case of two sets with large plane parallel faces 
of area A which are near to each other. In this case the associated mutual information has a leading term 
of the form 

I(A,B)~ k jp-z. (44) 
The coefficient k has been calculated for free fields in [1U|, 118] . In the extensive case, we have for odd D 

*« = (-D w ^ ( ^" v[5 L a) . («) 

and for D even 

7(J>/2-2)! 

KE 2 D / 2 - 1 (D-3)\l\(S D - 2 )' { ' 
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On the other hand, for the holographic entropy we have [6] 



D - D ~ X 



R D-1 2D-VD-D/2 /T^^y 

k »=^F^r{^ ] ^ (47> 

Then, we can compare the ratios of universal coefficients for the holographic and extensive functions 
in even dimensions 



D-l 

S^kE _ 1 / 1 ^(D-l)^ ' 

k H S E 2 D I 2 ^ D - 1 )/ 2 r ( D 



(48) 



Since this is different from one, it confirms the non-extensive character of the holographic function in this 
even dimensional case. 

4.6 Long distance limit 

Let us now look at the long distance limit of the mutual information. From (I15p it follows that for A and 
B contained in the same hyperplane, and in the conformal case, the extensive function is 

I (A, B) ~ Vol(A) Vol(fl) . (49) 

Note that this expression is invariant under independent rotations of A and B, and also independent on 
the direction of the vector of separation. 

Here we compare this behavior with the case of free massless scalars and fermions. These can be 
treated using the expression of the local density matrix in terms of correlators |28| . For a free fermion we 
have [23] 

2/3 



1/2 

where 



son- 



{R v ((3) - R v {-0)) 



(50) 



Rv(P)=[C v + p-h (51) 



2 

-1 

is the resolvent of the field correlator Cy restricted to V, 

C v (x, y) =< 0\*(x)**(v)\0 >= b(x -y) + tc^ Q \ X ~ V J D , (52) 

2 \x — y\ u 

where c is a constant depending on the dimension D and 7^ are the Dirac matrices. 

Considering V = A U B with a large separation vector L between A and B, the resolvent can be 
expanded perturbatively, 

Rvifi) = R V (P) ~ R%r(P)CiRviP) + R^m^RUWiRUP) ~ - • (53) 
Here Ry(f3) is the unperturbed resolvent, 

a$rW=( RA P (54) 



RbW) 

and C\ is the field correlator evaluated for the separation vector L 



The kernel 1a,b(x, y) = 1 for any x € A and y E B. 



. ( 1 A>B . 

c i = lc —nr- _i_ . n • (55) 
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The second term in (I53D does not contribute since it has zero trace. The third term is proportional to 
the squared of the field correlator and leads to 



I(A, B) ~ j^-q J d(3 ((3 - 1/2) [tr (r a V3)R%(P) + R B (0)R 2 A ((3)) - (0 - -(3)\ , (56) 

where V means the set V after an inversion of coordinates followed by a reflection in the plane perpendic- 
ular to L. The bar over the resolvent and the square of the resolvent means sum over the spatial variables, 

O x = fx dx fx d V °( x ' V)- 

Thus, the power I (A, B) ~ L -2 ^ -1 ) is the same as in the extensive case, but the coefficient may not 
be invariant under separate or joint rotations of A and B (it may depend on the direction of L). Also, 
the structure of this coefficient would in general be different to the product of the volumes. However, we 
have not been able to confirm this numerically in D = 3. 

A similar calculation for scalars shows that the mutual information also decays with the field correlator 
squared. For massless fields this gives I(A,B) ~ L _2 (- D_2 ) ) which is not consistent with extensivity. We 
can obtain this same result with a different approach. The mutual information is an upper bound on 
correlation [29] 

1 ((OaOb) - (Oa) (Ob)) 2 

2 I|Oa|| 2 ||Ob| 

where Oa and Ob are any bounded operators in the local algebras corresponding to A and B, with 
norms ||Oa|| and \\Ob\\ respectively. We cannot apply this relation directly to the smeared scalar field 
operator, j dx a(x)(j)(x), with a(x) an smooth function of compact support, since this is not bounded. 
However, using the inequality (|57|) with the unitary operators Oa = exp (i j dx a{x)(j){x)) and Ob = 
ex.p(i f dy P(y)4>(x)), and a(x), j3(x) vanishing outside A, B respectively, we obtain, 

I(A, B) > cons^(O)0(L)^ 2 ~ L" 2 ^- 2 ). (58) 

The same argument shows that for any hypothetical extensive theory there can not be correlators for 
localized bounded operators falling at a slower rate than L -2 ^ -1 ). 



77 A R\ ^ 1 W^A^B ~ \^A / R « 

I{A, B)>- — 2n ^ n2 , (57) 



5 Concluding remarks 

In conclusion, we have disproved extensivity of the mutual information for the Ryu-Takayanagi ansatz, 
as well as for the entropy functions for free scalar and Dirac fields (excepting the D = 2 massless Dirac 
case which is extensive). 

The deviations from extensivity, as measured by I(A,B,C), can have any sign. As far as we have 
checked, these tend to be relatively small (few percent) for the Dirac field, while they are much greater for 
the scalar field. The fact that different spins and sets configurations lead to different signs for I (A, B, C) 
gives some hope that a cancellation may occur for a particular interacting "extensive" theory. 

The entanglement entropy offers a non standard view of some aspects of QFT, and one can even 
hope that it might provide an alternative setting where QFT may be expressed in a geometrical way. An 
early attempt pursued this idea from an axiomatic point of view |llj . looking for the entropy functions 
S(V) which are subjected to the very general conditions of relativistic invariance, causality, positivity 
and strong subadditivity. The result can be reworded in the following terms. One either has divergent 
entanglement entropy, or the mutual information is zero. This last option is clearly inappropriate, since 
it implies the vanishing of all the correlators. 

Facing the fact that the entropy must diverge in a relativistic theory, one is lead to consider its universal 
part, which is encoded in the mutual information I(A,B). A list of requirements on this function is as 
follows: 

a) It is a symmetric (I(A,B) = I(B,A)) function on the pairs of spatially separated causally complete 
sets in Minkowki space-time 



12 



b) Poincare invariance 

c) Positivity 

d) It is monotonically increasing under inclusion 

e) I(A, B + x) goes to zero for large translations \x\ — > oo 

f) I(A, B, C) = I (A, B) + I(A, C) - I (A, B U C) is invariant under permutation of A, B and C. 

g) For two identical D — 1 dimensional cubes A and B (or equivalently their causal completion) with 
parallel faces, separated by a distance L, I(A,B) diverges as 

L -(D-2) ( Qr i og (£) for D = 2) in the short L 



We have added the condition (g) in order to allow only for solutions representing QFT with a finite 
number of fields. If this is not imposed one can construct solutions by dimensional reduction from higher 
dimension, as was done in Section 4. These types of functions have a more singular short distance behavior 
than the one given by (g) 0. 

Of course, any QFT gives a solution of the above (a-g) items. An interesting question is the converse. 
The expression (I27D with a positive, decreasing C(l), gives the first known explicit solution to this set 
of constraints. Another meaningful question is how to complete the holographic ansatz in order to fulfill 
these requirements. 
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